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ABSTRACT 

We present a proper motion mini-survey of 35 fields in the vicinity of Baade window, 
{l,b) = (1°,— 4°), sampling roughly a 5 x 2.5 deg region of the Galactic bar. Our 
second epoch observations collected with the ACS/HRC instrument on board the 
Hubble Space Telescope were combined with the archival WFPC2/PC images. The 
resulting time baselines are in the range of 4-8 years. Precise proper motions of 15,863 
stars were determined in the reference frame defined by the mean motion of stars with 
magnitudes between If814W = 16-5 — 21.5 along the line of sight. We clearly detect 
small gradients in proper motion dispersions (ci, crb) ^ (3.0, 2.5) mas yr^^, and in the 
amount of anisotropy (cri/(Tb ~ 1.2). Both the longitude dispersion ci and its ratio to 
the vertical motion Cb increase toward the Galactic plane. The decline of the anisotropy 
ratio cri/(Tb toward the minor axis of the bulge is mostly due to increasing ah- We also 
find, for the first time, a significant negative covariance term in the transverse velocity 
field (T\h / {(J\<7h) — —0.10. Our results extend by a factor of ~15 the number of the 
Galactic bar fields with good proper motion dispersions. 

Key words: galactic dynamics - Stars: proper motion, dispersion 
gravitational lensing - Galaxy: bar, bulge, disc 



1 INTRODUCTION 

The Milky Way appears to be a typical spiral galaxy with a 
disk and bulge. While our unique inside view of the Galaxy 
helps to understand the galactic structure in general, it also 
makes it more difficult to identify structures such as bars. 
The case fo r existence of a bar at the Galactic centre - first 
proposed bv lde Vaucouleura l|l96^ - is easier to make know- 
ing that bars are common in external galaxies. 

There is now conclusive evidence that the Galac- 
tic bulge is of a barred type. Th e longitud e asymme- 
try of the COBE photometric maps JBlitz fc Spergeli,199ll 



* Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc. under NASA 
contract NAS 5-26555. 
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iDwek et al.l 1995|), high optical depths t o gravitational mi 
crolensi ng Zhao^gerge^&^£hj_^95|)_j asymmetric star 



counts iStanek et al.l Il994l : jBaMsiau^^j^Gilmord 12005 
non-circular gas kinematics ide Vaucou lcurs" 1964), and tri- 
axial i ty of the stellar velocity f ield iZhao. Spcrgcl fc RichI 
19941: IZhao. Rich fc Bieiiolll996^ have aU been interpreted 
as signatures of the Galactic bar. Unfortunately, the size 
and preci se orientation o f the b ar are still being debated. 
Recently iBeniamin et alJ ^2005^ found that the infra-red 
star counts collected by the Spitzer Space Telescope are 
best explained assuming a bar with a half-length 4.4 ± 0.5 
kpc placed at a ~ 44° angle to the Sun-Galactic center 
line. Most previous studi es prefer a bar at ~ 20° extend- 
ing out to ~3.5 kpc (e.g. lGerhardll2001f) . Such apparently 
conflicting evidence may be an indication that the inner 
Gala xy hides even more complicated structures. A secondary 
bar ( Alard '2001'; 'Babusiaux fc Gilmor'3 l2005l) and a ring 
(iSevenster fc Kalnajs 2001) have been suggested, since these 
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features are also evide nt in many other spiral galaxies (e.g. 
lErwin fc Sparkdll999h . 

iBinnevI 112005^ discussed the progress of the dynam- 
ical modeling techniques in the context of major ob- 
servational advances expected from a future space mis- 
sion GAIA. The three approaches to constructing a self- 
consisten t dynamical Galaxy model a re the Schwarzschil d 
method JZhao. Soergel fc RichI >1994': Hafncr et alj l20o3) . 
the torus modelling method (see Bimicv 20051 for details) 
and A''-body simulations with partic le weights determined b y 
the "made-to-measure" algorithm JSver fc Tremainelll99£i) . 
The first Gal actic bar mode l empl oying the latter method 
was built by iBissantz et al.l i2004ft . Neither of these tech- 
niques can fully address the structure of the inner Galaxy 
without constraints on stellar kinematics. The refinement of 
the models is limited largely by the scarcity of good proper 
motion and radial velocity measurements. Bissantz ct al. 
l|2004) . for example, compared kinematic predictions of their 
model with the data for just two lines of sight. A handful of 
samples published since t he pioneering photographic work of 
ISpaenhauer et al.l l|1992h is not enough to remove the non- 
uniqueness of the model parameters. 

In a study based on two lines of sight iKuiiken fc RichI 
l|2002l) have demonstrated that high quality relative proper 
motions can be obtained with a relatively modest investment 
of time using the Hubble Space Telescope (HST). At the 
resolution of the Wide Field Planetary Camera 2 (WFPC2) 
instrument the required time baseline is only a few years. 
The HST archive contains a number of images suitable as 
the first epoch data, so the tedious part of accumulating 
the baseline can be avoided entirely. Most of these fields are 
centred around microle nsing events discovere d by the MA- 
CHO collaboration (e;gJPo£Owsk^^L|2005|) . Using a sim- 
ilar concept to that of lKuiiken fc Rlchl l|2002h . we carried out 
a mini-survey of proper motions in 35 of the available MA- 
CHO fields to study the kinematics of microlensed sources 
and of the general stellar populations. Here we present the 
measurement techniques and results for the general stellar 
population in these 35 lines of sight. 



2 HST IMAGES AND DATA REDUCTION 

The log of the HST observations used in our proper mo- 
tion mini-survey is given in Tabled The first epoch images 
(selected from the HST archive ) were all taken with the 
WFPC2/PC camera, and cover the time interval 1996-2000. 
The second epoch data come from our own SNAP program 
(cycle 13; proposal ID 10198) and were collected in 2004 
and 2005 using the High Resolution Channel (HRC) of the 
Advanced Camera for Surveys (ACS). Our SNAP survey 
was optimized toward high execution rates and, therefore, 
we only requested F814W observations to keep the required 
target visibility as low as possible. Both PC and HRC detec- 
tors cover a similar field of view (25" x 29" and 35" x 35", 
respectively) and have comparable pixel scales (27 versus 
45.5 mas). There were no restrictions on the telescope roll 
angle during ACS observing. While the latter relaxed con- 
dition decreased the number of possible proper motion de- 



terminations, it greatly improved scheduling opportunities. 
Most of the subsequent analysis for each of the 35 fields 
is based on a pair of good quality F814W (7-band) images 
constructed by stacking all suitable data for a given epoch. 
In some cases, the first epoch data included F555W {V- 
band) images that allowed us to construct color-magnitude 
diagrams (CMD s). We also re-analyzed the two fields previ- 
ously studied bv lKuiiken fc R ich ( 2002), increasing to 37 the 
total number of the Galactic bulge fields considered here. 



2.1 Image reductions 

The basic reductions of the ACS images, i.e. de-biasing, 
dark frame subtraction, flat-fielding, and cosmic-ray re- 
moval, were performed on-the-fly by the standard HST data 
processing pipeline. The pipeline also takes care of dither- 
ing, cosmic-ray splits and geometric corrections using the 
Multidrizzle software ( Kookcmocr ot al. 2002) , which in turn 
uses the Drizzle routines (Fruchter fc Hook 2002). Our ACS 
observations employed a generous 4-point dithering pattern 
combined with a 2- way cosmic- ray split, providing the fi- 
nal drizzled images with high S/N ratio, excellent dynamic 
range and highly reliable cosmic-ray rejection. In case of the 
first-epoch WFPC2 images we utilized the standard HST 
data products for individual exposures, and then used the 
drizzle task of the IRAF package to correct the geometric 
distortions. For cosmic-ray cleaning, registering and com- 
bining these corrected images we developed dedicated IRAF 
scripts. The quality of our final cross-instrument astrometry 
is limited by the larger pixel size, as well as the lower S/N 
ratio and number of the individual first epoch PC frames 
available for stacking by comparison to the ACS data (see 
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2.2 Object catalogs and PSF fitting 

The instrumental positions and magnitudes of the field 
objects were measured using the IRAF task starfind, 
an improved version of daofind that fits Gaussian pro- 
files to stellar images. The combined images from both 
WFPC2/PC and ACS/HRC detectors have a well sam- 
pled point-spread-function (PSF) with the fuU-width-at- 
half-maximum (FWHM) of stellar images, correspondingly, 
2.4 and 2.8 pixels. Our PSF fits were restricted within the 
area of the Airy disc (3.0 and 2.0 pixel radius for PC and 
HRC data respectively), where the point source fiux is well 
approximated by a Gaussian model. Outside the Airy disc 
the PSFs show a variety of shapes, including rings, possible 
diffraction spikes and bright spots in case of high S/N ob- 
jects. These features can mimic stars and need to be care- 
fully considered during object cross-identification. We im- 
posed a minimum separation of 3 x FWHM between any two 
sources detected in the same image to ensure that there are 
essentially no spurious objects in the final source lists. The 
loss of number statistics due to the accidental rejection of 
the actual stars in tight groups is insignificant. In fact, the 
centroid measurements for objects in the wings of other stars 
are notoriously unreliable and best avoided. The minimum 
separation cut also helps in cross identification of objects 
between the two epochs (§|21l, since the expected intrinsic 
object shifts may reach ~2 PC pixels. 



© 2006 RAS, MNRAS 000, 000-000 



Table 1. Log of the HST observations. 



First epoch 



Second epoch 



MACHO field 


RA 


Dec 


Year 


F814W exp. 


F555W exp. 


Prop. ID 


Year 


F814W exp. 


108-C 


18:00:01.276 


-28:27:41.23 


1996.82 


6 X 260 s 




6756 


2005.16 


4 X 260 s 


119-C 


18:03:03.010 


-30:09:56.50 


1996.82 


6 X 260 s 




6756 


2005.15 


4 X 260 s 


119-D 


18:04:24.825 


-30:05:58.94 


1996.82 


6 X 260 s 




6756 


2004.78 


4 X 260 s 


120-A 


18:07:26.441 


-29:39:34.22 


1996.82 


6 X 260 s 




6756 


2005.15 


4 X 260 s 


167-A 


18:13:32.154 


-26:31:10.33 


1996.82 


6 X 260 s 




6756 


2005.16 


4 X 260 s 


101-C 


18:07:32.649 


-27:31:35.60 


1997.47 


6 X 260 s 




6756 


2005.15 


4 X 260 s 


95-BLG-ll .. 


18:04:37.239 


-30:12:11.45 


1996.67 


6 X 260 s 




6756 


2005.16 


4 X 260 s 


96-BLG-17 .. 


18:06:09.107 


-27:53:38.59 


1996.81 


6 X 260 s 




6756 


2005.17 


4 X 260 s 


119-A 


18:03:35.789 


-29:42:01.26 


1996.68 


6 X 160 s 


2 X 400 s 


6756 


2005.14 


4 X 160 s 


95-BLG-7 ... 


18:13:29.298 


-26:13:58.12 


1998.84 


2 X 40 s 


3 X 40 s 


7431 


2005.46 


4 X 40 s 


95-BLG-lO .. 


17:58:16.011 


-29:32:10.86 


1997.82 


2 X 40 s 


3 X 40 s 


7431 


2004.66 


4 X 40 s 


95-BLG-13 .. 


18:08:47.038 


-27:40:47.25 


1999.45 


2 X 40 s 


3 X 40 s 


7431 


2005.12 


4 X 40 s 


95-BLG-14 .. 


18:01:26.308 


-28:31:14.03 


2000.45 


2 X 40 s 


3 X 40 s 


7431 


2005.40 


4 X 40 s 


95-BLG-19 .. 


18:11:32.487 


-27:45:26.99 


1998.49 


2 X 40 s 


3 X 40 s 


7431 


2005.44 


4 X 40 s 


97-BLG-18 .. 


18:03:15.254 


-28:00:14.06 


1998.59 


2 X 40 s 


3 X 40 s 


7431 


2005.31 


4 X 40 s 


104-C 


18:03:34.050 


-28:00:18.94 


1998.73 


2 X 40 s 


3 X 40 s 


7431 


2005.43 


4 X 40 s 


104-D 


18:03:29.024 


-28:00:30.99 


1998.80 


2 X 40 s 


3 X 40 s 


7431 


2005.45 


4 X 40 s 


108-A 


18:00:25.866 


-28:02:35.24 


1998.80 


2 X 40 s 


3 X 40 s 


7431 


2005.16 


4 X 40 s 


128-B 


18:07:18.624 


-28:59:29.83 


1998.49 


2 X 30 s 


3 X 40 s 


7431 


2005.37 


4 X 30 s 


104-B 


18:03:09.046 


-28:01:45.25 


1999.45 


2 X 40 s 


3 X 40 s 


7431 


2005.39 


4 X 40 s 


128-A 


18:06:57.621 


-29:00:55.15 


1999.33 


2 X 40 s 


3 X 40 s 


7431 


2005.49 


4 X 40 s 


94-BLG-3 ... 


17:58:25.300 


-29:47:59.50 


1997.82 


2 X 40 s 


3 X 40 s 


7431 


2005.48 


4 X 40 s 


94-BLG-4 ... 


17:58:36.766 


-30:02:19.27 


1997.82 


2 X 40 s 


3 X 40 s 


7431 


2005.16 


4 X 40 s 


95-BLG-36 .. 


18:07:20.775 


-27:24:09.69 


1998.80 


2 X 40 s 


3 X 40 s 


7431 


2005.39 


4 X 40 s 


95-BLG-37 .. 


18:04:34.452 


-28:25:33.46 


1999.43 


2 X 40 s 


3 X 40 s 


7431 


2004.67 


4 X 40 s 


95-BLG-38 .. 


17:59:41.851 


-28:12:10.31 


1998.81 


2 X 40 s 


3 X 40 s 


7431 


2005.33 


4 X 40 s 


95-BLG-41 .. 


18:02:06.332 


-28:50:45.26 


1999.46 


2 X 40 s 


3 X 40 s 


7431 


2005.44 


4 X 40 s 


96-BLG-14 .. 


18:05:15.421 


-27:58:25.01 


1997.83 


2 X 40 s 


3 X 40 s 


7431 


2004.67 


4 X 40 s 


96-BLG-4 ... 


18:06:11.954 


-28:16:52.77 


1998.79 


2 X 26 s 


3 X 40 s 


7431 


2004.82 


4 X 26 s 


97-BLG-38 .. 


18:04:06.083 


-27:48:26.25 


1998.51 


2 X 40 s 


3 X 40 s 


7431 


2004.63 


4 X 40 s 


97-BLG-24 .. 


18:04:20.253 


-27:24:45.28 


1998.35 


2 X 40 s 


3 X 40 s 


7431 


2005.49 


4 X 40 s 


96-BLG-5 ... 


18:05:02.497 


-27:42:17.23 


1999.45 


4 X 160 s 


2 X 400 s 


8490 


2005.12 


4 X 160 s 


98-BLG-6 ... 


17:57:32.812 


-28:42:45.41 


2000.48 


2 X 100 s 


2 X 260 s 


8654 


2004.73 


4 X 100 s 


97-BLG-41 .. 


17:56:20.691 


-28:47:41.97 


2000.47 


4 X 100 s 


4 X 160 s 


8654 


2004.62 


4 X 100 s 


99-BLG-22 .. 


18:05:05.281 


-28:34:41.69 


2001.77 


4 X 400 s 


4 X 400 s 


9307 


2005.16 


4 X 400 s 



The final object catalog s were conv erted to the 
VEGA magnitude system (Girardi et alj2002l and references 
therein) and the astrometric transformations to the Galactic 
(1, b) coordinates were established using the World Coordi- 
nate System (WCS) headers of the ACS images. Our esti- 
mated S/N ratios for object fluxes are based on propagated 
errors in pixel counts that account for photon statistics. 



3 ESTIMATING TRANSVERSE MOTIONS OF 
THE GALACTIC BULGE STARS 

Absolute astrometry is difficult in the crowded Galactic 
bulge fields. Until we can establish a sample of extragalactic 
objects (e.g. spectros copically confirme d QSOs in the cata- 
logue of candidates bv lSumi et al.l2005h shining through the 
low extinction windows, the only readily available reference 
velocity in the Galactic bulge is the mean velocity of stars 
along the line of sight. Notice, however, that the second order 
moments of proper motions are unaffected by the choice of 
our reference frame. In this analysis we use the magnitude- 
and distance-selected samples to investigate the spatial de- 



pendence of the covariance matrix of the transverse velocity 
field across the Galactic bar. 



3.1 Relative proper motions and their dispersions 

Having measured the instrumental positions of stars on 
both the first and the second epoch images, we tied the 
WFPC2/PC positions to the ACS/HRC pixel grid. The ob- 
ject shifts Al and Ab in the Galactic coordinates between 
the two epochs could then be calculated using the WCS in- 
formation from the ACS headers. We cross-correlated the 
positions of a few hundred stars in the magnitude range 
16.5 < /f814w < 21.5 to obtain the coordinate transforma- 
tion between two pixel grids, which is approximated by a 
third order polynomial. Stars brighter than 7f8i4w = 16.5 
were often saturated while those with 7f8i4w > 21.5 were 
too faint to have useful S/N ratios, particularly for the fields 
with short exposures (Table . Our procedure for cross- 
identifying stars starts from matching the first 20 objects 
(out of ~50 brightest sta rs with 17 < Jf814w < 18 ) us- 
ing the triangle algorithm (lGrothll986t IWozniakl2000ll . The 
initial low-order fit is then iteratively refined. A star with a 
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transverse velocity of 100 km/s at the distance of 8 kpc will 
move by 26.4 mas, or roughly one ACS/HRC pixel, assum- 
ing a 10 year baseline. Accordingly, we adopted a tolerance 
radius of 100 mas for the final matching. 

After geometrically aligning and transforming object 
positions to the Galactic (1, b) coordinates, we folded the 
data with the time baseline and estimated all components of 
the transverse velocity tensor, i.e. dispersions cti, ab and the 
normalized covariance Cib = crib/(o"iab). The sample of stars 
used to trace the kinematic parameters of the Galactic bulge 
was limited to the magnitude range 18.0 < /f814w < 21.5, 
i.e. dominated by the bulge main sequence population near 
the turn-off point. This puts all lines of sight (with data sets 
of the varying depth and dynamic range) on a more com- 
mon f ooting. However, as already noted bv lKuiiken fc RichI 
l)2002r . the results are insensitive to the details of the mag- 
nitude cuts. 



3.2 Astrometric errors 

The la centroid errors from PSF fitting (per coordinate) 
can be estimated from the S/N ratio (SNR): 



: 7 X 



FWHM 
SNR ' 



(1) 



where 7 = 0.6 fo r a Gaussian PSF mod el and the FWHM is 
in pixels (see e.g. lKuiiken fc Richll2002h . We tested this pre- 
scription by stacking independent subsets of images taken at 
a single epoch. Similarly to Kuiikcn & Rich (2002), we find 
that Equation is an excellent representation of the actual 
astrometric uncertainties in our data, with the exception 
of the brightest stars, for which a constant systematic con- 
tribution of 0.025 pixel is required. Consequently, we used 
Equation Q with the systematic term added in quadrature 
to estimate the astrometric errors and their contribution to 
the apparent proper motion dispersions. The formulas for 
estimating cti, (Jb and their erro rs corrected for the mea - 
surement variance can be found in S oaenhaxier et al.l l|1992ll . 
Throughout this paper we use bootstrapped uncertainties of 
the sample statistics (from 1000 trials) that turned out to be 
slightly more conservative than analytical formulas. The es- 
timated intrinsic dispersions reported in §^are 5-10% lower 
compared to the raw values. The cross term Cib need not be 
corrected, as long as the errors in 1 and b are uncorrelated. 
None of our conclusions depend on the precise value or even 
the presence of this correction. 

The limiting S/N ratio for a useful detection in our anal- 
ysis is about 10 and corresponds to a /f8i4w — 21.5 mag 
star in the combined image of two 40-second WFPC2/PC 
exposures. The same star will be detected at S/N ~ 20 in 
the lowest quality ACS stack (four 40-second frames). The 
shortest time baseline in our data is 3.388 years, and the 
typical la astrometric uncertainties for a 21.5 mag star are 
~7.2 and ~2.1 mas in the first and the second epoch images 
respectively. In this worst case scenario, the proper motion 
can be measured to an accuracy 2.5 mas yr~^. The images 
for the first eight fields in Tabled have relatively long ex- 
posure times, so the resulting proper motion errors are only 
~0.1 mas yr~^ for bright stars and ~0.3 mas yr~^ for the 
faintest stars in those samples, with a systematic error of 
0.025 pixel (c.f. discussion following equation Q). 



4 RESULTS 

The results are given in Table |5| and plotted in Figs. and 
121 Proper motions for individual stars are available online'^. 
After presenting our measurements we check for consistency 
with two other published data sets (S I4.3II . A more detailed 
discus sion and comparison to the results of lKuiiken fc RichI 
J2002D follows in §|11 



4.1 Proper motion dispersions 

The spatial dependence of proper motion dispersions ai, ah 
is shown in Fig.Q Recall that at a distance of 8 kpc, a veloc- 
ity of 100 kms~^ implies a proper motion of 2.64 mas yr~^. 
The most visible trends are in cri(b) and crb(l) that tend to 
increase closer to the Galactic plane and the Galactic centre. 
Both gradients are weak, but clearly present. From a simple 
straight line fit we find: 



ai = 0.16 ± 0.04 X b + 3.38 ± 0.13, 



fT-b 



-0.09 ±0.02 X 1 + 2.62 ±0.06. 



(2) 
(3) 



Consequently, ah increases from 2.1 to 2.6 mas yr~^, or by 
about 20%, as the longitude 1 varies from 5.5 to 0.5 deg. 
Similarly, ah changes from 2.6 to 3.2 mas yr~^ between b = 
—4.5 and —2 deg. It is intriguing that the last data point 
around fe ~ — 2 deg has the lowest dispersion (ai) measured 
for all fields (see the top right panel in Fig. 0, but the 
value is still marginally consistent with the observed scatter. 
There is no other indication of the intrinsic variations on the 
field-to-field scale. The distributions of cri(l) and (Jb(b) are 
fiat within the scatter from random errors and projection 
effects in the presence of gradients. 



4.2 Anisotropy and covariance 

In Fig. 121 we plot the ratio and the correlation coefficent 
(covariance of the velocity field) of cti and ah as a func- 
tion of location in the bulge. There is a significant level 
of anisotropy, i.e. ai/ah > 1, throughout the covered area. 
Moreover, the velocity distribution shows a tendency to be- 
come more isotropic for lines of sight approaching the Galac- 
tic centre at a roughly fixed latitude b. This is a reflection 
of the increase in ah with approximately constant ai (S I4.1II . 
The trend of more anisotropy toward the Galactic plane is 
also driven primarily by one of the dispersions (cri), but it is 
more difficult to see. Part of the reason for this is the narrow 
range of b covered by the data. The formal fits give: 



ai/ah = 0.05 ± 0.01 x 1 ± 1.08 ± 0.03, 
ai/ah = 0.03 ± 0.03 x b ± 1.27 ± 0.08. 



(4) 
(5) 



The estimates of the covariance term from Table [51 (plot- 
ted in Fig. 121 are all negative and scatter uniformly in the 
range —0.20 < Cib < —0.02. This indicates that in our 
Galactic bulge fields the stellar motions in directions par- 
allel and perpendicular to the plane are significantly anti- 
correlated. An observational bias that would account for 
this anti-correlation has to operate in a similar way over 
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Figure 1. Spatial dependence of proper motion dispersions cti and cri, in Galactic coordinates for our turn-off point dominated sample 
in the Galactic bulge ("Table l2l and Sl4l. The two open circles are for the Baade window and Sagittarius-I fields from Kuijken &; Rich 
(2002). The lines show linear regressions (solid) and weighted means (dashed) of the data. For the top right panel, the rightmost data 
point was not used in the fit. 







0.1 


-0.2 
-0.3: 



1 1 1 1 1 


1 1 1 1 1 

1 1 1 1 1 


1 1 1 1 1 



4 3 2 

1 [deg] 



1 








b [deg] 

Figure 2. Similar to Fig. [Tl but for the anisotropy ratio cri/cr^, and covariance term Cib = o-i\^ / {aio-^^) . 
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Table 2. Results of our proper motion mini-survey. The dispersions cri,ai, and the dimensionless 
correlation coefficient Cib were measured for 35 lines of sight in the Galactic bulge (1, b). The time 
baseline At and the number of stars Avatars used to estimate the kinematics are also given. 



Field name 



1 
[deg] 



b 
[deg] 



mas yr" 



Cb 
[ mas yr^-"-] 



Cn 



At 

[yr] 



iVsl 



101-C 


3.65 


-3.47 


2.85 ±0.09 


2.45 ± 0.08 


-0.15 ±0.05 


7.683 


445 


104-B 


2.73 


-2.87 


2.97 ±0.10 


2.50 ±0.10 


-0.05 ±0.05 


5.941 


407 


104-C 


2.80 


-2.93 


2.74 ± 0.09 


2.51 ±0.10 


-0.15 ±0.04 


6.706 


482 


104-D 


2.79 


-2.92 


2.84 ±0.10 


2.36 ±0.10 


-0.10 ±0.05 


6.649 


437 


108-A 


2.42 


-2.35 


2.90 ±0.12 


2.32 ±0.12 


-0.08 ±0.06 


6.360 


396 


108-C 


2.02 


-2.48 


3.15 ±0.10 


2.52 ±0.07 


-0.09 ± 0.04 


8.345 


615 


119-A 


1.32 


-3.77 


2.89 ±0.10 


2.44 ± 0.08 


-0.14 ±0.04 


8.458 


471 


119-C 


0.85 


-3.89 


2.79 ±0.10 


2.65 ±0.08 


-0.14 ±0.04 


8.339 


459 


119-D 


1.06 


-4.12 


2.75 ±0.10 


2.56 ±0.09 


-0.05 ±0.06 


7.962 


420 


120-A 


1.76 


-4.48 


2.75 ±0.09 


2.52 ±0.09 


-0.04 ±0.05 


8.339 


397 


128-A 


2.28 


-4.08 


2.63 ±0.11 


2.33 ±0.12 


-0.12 ±0.05 


6.165 


357 


128-B 


2.33 


-4.13 


2.70 ±0.12 


2.29 ±0.13 


-0.13 ±0.06 


6.881 


338 


167-A 


5.17 


-4.16 


2.75 ±0.11 


2.36 ± 0.09 


-0.18 ±0.05 


8.345 


317 


94-BLG-3 .. 


0.68 


-2.84 


2.84 ±0.10 


2.58 ±0.10 


-0.12 ±0.05 


7.654 


496 


94-BLG-4 .. 


0.49 


-3.00 


2.58 ±0.11 


2.46 ± 0.09 


-0.03 ± 0.04 


7.341 


413 


95-BLG-lO . 


0.89 


-2.68 


3.07 ±0.10 


2.41 ±0.09 


-0.12 ±0.04 


6.840 


487 


95-BLG-ll . 


0.99 


-4.21 


2.82 ± 0.09 


2.62 ±0.09 


-0.14 ±0.04 


8.493 


443 


95-BLG-13 . 


3.64 


-3.78 


2.61 ±0.13 


2.31 ±0.12 


-0.14 ±0.05 


5.672 


309 


95-BLG-14 . 


2.12 


-2.78 


2.95 ±0.13 


2.50 ±0.11 


-0.12 ±0.05 


4.950 


463 


95-BLG-19 . 


3.87 


-4.36 


2.61 ±0.11 


2.17±0.10 


-0.13 ±0.06 


6.952 


300 


95-BLG-36 . 


3.73 


-3.37 


2.75 ±0.12 


2.11±0.11 


-0.10 ±0.05 


6.587 


376 


95-BLG-37 . 


2.54 


-3.33 


2.72 ±0.12 


2.44 ±0.12 


-0.04 ±0.05 


5.238 


442 


95-BLG-38 . 


2.20 


-2.29 


2.87 ±0.12 


2.46 ±0.10 


-0.05 ± 0.04 


6.526 


474 


95-BLG-41 . 


1.91 


-3.07 


2.79 ±0.10 


2.34 ±0.10 


-0.04 ±0.05 


5.980 


450 


95-BLG-7 .. 


5.42 


-4.01 


2.86 ±0.14 


1.88 ±0.11 


-0.20 ±0.07 


6.616 


265 


96-BLG-14 . 


3.01 


-3.24 


2.71 ±0.12 


2.40 ±0.12 


-0.17 ±0.05 


6.833 


373 


96-BLG-17 . 


3.17 


-3.38 


3.07 ±0.10 


2.55 ±0.09 


-0.16 ±0.04 


8.364 


557 


96-BLG-4 .. 


2.84 


-3.57 


2.68 ±0.14 


2.26 ±0.14 


-0.04 ±0.06 


6.027 


329 


96-BLG-5 .. 


3.22 


-3.07 


3.17 ±0.10 


2.39 ± 0.08 


-0.13 ±0.05 


5.670 


535 


97-BLG-18 . 


2.77 


-2.87 


2.99 ±0.10 


2.38 ±0.10 


-0.12 ±0.04 


6.713 


433 


97-BLG-24 . 


3.40 


-2.79 


3.00 ±0.11 


2.39 ±0.10 


-0.10 ±0.05 


7.115 


398 


97-BLG-38 . 


3.03 


-2.94 


2.95 ±0.12 


2.21 ±0.10 


-0.06 ±0.05 


6.118 


395 


97-BLG-41 . 


1.32 


-1.95 


2.58 ±0.07 


2.13 ±0.07 


-0.09 ±0.04 


5.145 


612 


98-BLG-6 .. 


1.53 


-2.13 


3.26 ±0.10 


2.79 ±0.12 


-0.07 ±0.05 


4.252 


670 


99-BLG-22 . 


2.46 


-3.50 


3.11 ±0.10 


2.60 ±0.09 


-0.17 ±0.04 


3.388 


493 


KR-BW^ .. 


1.14 


-3.77 


2.87 ±0.08 


2.59 ±0.08 


-0.07 ±0.03 


6.048 


694 


KR-SgrI ... 


1.26 


-2.66 


3.07 ±0.08 


2.73 ±0.07 


-0.09 ±0.04 


5.960 


752 



'^ KR -BW and KR-SgrI are the Baade window and the Sagittarius-I field from iKuiiken fc RichI 
<2002l) 



a large range of instrumental settings. For example, a seri- 
ous concern is a presence of preferred telescope orientations. 
Indeed, for about half of our fields the relative roll angle be- 
tween the two compared observations falls in a narrow range 
of 25 deg. The other half, however, is spread over all possi- 
ble orientations and still shows about the same covariance. 
The skewness of the ACS focal plane cannot be the cause of 
the observed correlatio n, because the measurements in both 
iKuiiken fc RichI ^QQ^ fields use only WFPC2/PC data and 
yet they perfectly agree with the rest of the Cib values. We 
also investigated several other possibilities, we found no ex- 
planation for this result other than a true correlation be- 
tween /ii and ^h- Taking a field with a relatively low S/N 
ratio in our data and assuming perfectly correlated errors 
in /ii and ^h, the expected covariance is only Cib ~ 0.02. 



There is a slight hint in Fig. |5]that Cib may vary with the 
longitude, although this impression seems to rely on the two 
points farthest from the bulge minor axis (1~5.3 deg in the 
left panel of Fig. |2j that fall below the rest of the data. 



4.3 Comparisons with previous work 

4-3.1 OGLE-II proper motion cat 



used large number 
( U dalski et al. 199^ 



statistics of the 
to derive relative 



ISumi et al.l ||21 
OGLE-II database 
proper motions of >5xlO'' stars in the Galactic bulge re- 
gion from hundreds of observations covering a 4- year base- 
line. The OGLE-II catalogue is a valuable resource for kine- 
matic studies of bright stars like red clump giants that are 
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=1. 




^,,„,^^[mas/yr] 



[mas/yr] 



Figure 3. Comparison between our HST measurements and the ground based OGLE-II data for bright stars from the catalogue of Sumi 
et al. (2004). There are 77 stars covered by our observations that have catalogue errors 3 mas yr~^ or better in Sumi et al. Significant 
discrepancies (marked by alphabets) are caused by blending (c.f Fig. |1J. The solid lines indicate the two measurements (ground-based 
and from the HST) are equal. 



relatively free of the source confusion effects. However, at 
the 1.3" FWHM seeing of the ground based OGLE-II im- 
ages, a random red clump giant star still has a few per cent 
probability of being blended with another unre solved star. 
It is i nstructive to cross- validate the results of ISumi et alJ 
l|2004) and our high-resolution HST measurements against 
each other. 

Out of 35 program fields in Tabled 15 are covered by 
the OGLE-II proper motion catalogue. In our HST sample 
we found 77 stars for which the OGLE-II catalogue proper 
motion error is ^ 3 mas yr~^. The two data sets were com- 
pared star by star after adjusting for an arbitrary zero point 
of the proper motion scale. The results are plotted in Fig.|3] 
and show a good overall agreement of our measurements 
with Sumi et al. (2004). All significant outliers were labeled 
and checked for blending. Fig.^Jdemonstrates that virtually 
all these substantial discrepancies are linked to the presence 
of an unresolved companion within ~1" of the primary ob- 
ject. 



i.3.2 Kmjken & Rich (2002) 

Our approach to measuring the positions and proper mo- 
tions o f stars (SI2.2II is somew hat simpler than the method 
used bv lKuiiken fc Richlll2002h . The latter study utilized the 
images from the WFPC2/WF chips and had to accommo- 
date a strong under-sampling of the PSF. In contrast, our 
use of the WFPC2 data was limited to the critically sampled 
images from the PC detector. The second epoch ACS/HRC 
images have four times the PSF sampling of the WF images, 
so we could take advantage of the conventional PSF fitting 
techniques. 

Regardless, in order to eliminate the possibility of a hid- 
de n error we re-analysed the PC data in both fields studied 
bv lKuiiken fc RichI i2002r using our tools. Table 3 shows the 
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Figure 4. Cutout HST images for the outliers marked by al- 
phabets in Fig. 13 All significant outliers in Fig. Islcan be linked 
to source confusion and flux blending. The dark circles have a 
diameter of 1 arc second. 



results of this comparison. The agreement between the two 
sets of measurements is remarkably close despite significant 
differences in the sample size and the adopted selection cri- 
teria. This also confirms that our results are not significantly 
affected by several subtle instrumental effects that can po- 
tenti ally infiuence astrome tric work with the HST images 
(e.g. iKuiiken fc Richll2002l and references therein) . 
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Table 3. Proper motion dispersions from Kuijken & Rich (2002) compared with the results of our reanalysis of the 
same data 





1 


b 




This work 




Kuiiken & Rich ('2002") 




Field 


o-i 


Cb 


A^stars 


o-l 


o-b 


A^stars 


BW 
Sgr-I 


1.14 
1.27 


-3.77 
-2.66 


2.87 ± 0.08 
3.07 ± 0.08 


2.59 ± 0.08 
2.73 ± 0.07 


694 
752 


2.91 ± 0.06 
3.10 ± 0.06 


2.51 ± 0.05 
2.73 ± 0.05 


1076 
1388 
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Figure 5. Color-magnitude diagram (CMD; left) and relative proper motions (right) for stars in three nearly coincident stellar fields 
from Table 121 97-BLG-18, 104-C and 104-D. The red and blue stars above the turn-off point show the kinematics characteristic of the 
bulge and disk populations, respectively. 
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Figure 6. Histograms of relative proper motions of the red (solid line) and blue (dashed line) samples from Fig. |F] The blue disk stars 
"rotate in front" of the Galactic bulge parallel to the plane. 
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Figure 7. Average relative proper motions and dispersions of stars in the Baade Window in bins of M* = /f814W — 2 x (Vfsssw ~ ^F814W ) ; 
an approximate distance indicator. 



5 DISCUSSION 

5.1 Distance and population trends 

The study of lKuiiken fc Ricbl i2002r focused on cleaning the 
Galactic bulge population and removing the contamination 
by the bluer disk stars. Above the bulge turn-off point, the 
stellar colors alone are sufficient to separate the blue disk 
main sequence from the red giants, subgiants and clump 
giants. The size of our fields is generally too small to pro- 
vide useful statistics of bright stars above turn-off point, and 
good color information is only available for about 1/3 of the 
lines of sight. However, three of the fields in Tabled with 
useful colors (97-BLG-18, 104-C and 104-D) are close to each 
other and were combined in order to look for a kinematic 
distinction between the Galactic disk and bulge populations. 
Figs 13 and |S| show that in the frame of reference of a mean 
star (of any color), the longitude proper motions of the blue 
disk stars are biased toward positive values, while the red 
bulge stars tend to have more negative ni. The blue and 
red samples were selected, correspondingly, using conditions 
(V - I) < 1.65 and (V - I) > 1.7 . This effect was previ- 
ously observed bv lSumi et alj ^2005^ and bv lKuiiken fc RichI 
i2002l) in their two fields with multi-epoch WFPC2 data and 
there is little doubt that it is due to the blue disk stars "ro- 
tati ng in front" of the red bulge stars. 

iKuiiken fc RichI i2002l) also devised an approximate dis- 
tance measure: 



M = /f814W — 2(Vf555W — -^F814w) 



(6) 



chosen to remove the slope of the main sequence in the 
color- magnitude diagram. In Fig. |7|we present the average 
proper motions and their dispersions for the Baade window 
in bins of Af *. As expected, with an increasing depth along 
the line of sight, the kinematic signature gradually changes 
from that characteri stic of the disk stars, to the one typical 
for the bulge. In the lKuiiken fc RichI i2002|) data this trend 
continues to very faint stars that are likely on the far side 



of the bulge, and if so, it constitutes a "rotation curve" of 
the bulge. The colors for our fields are generally of lower 
S/N ratio or nonexistent, and do not allow to see this much 
detail. 



5.2 Stellar velocity ellipsoid of the Galactic bar 

A detailed modeling of the measurements in Table |2| is be- 
yond the scope of this paper. Here we only comment on 
pos sible directions and n e w poss ibilities. 

IZhao. Rich fc Bielld il996l) interpreted the bulge 
anisotropy in terms of the rotation support of the Galac- 
tic bulge and related the ratio ai/ah directly to the level of 
flattening of the light density distribution. They also con- 
cluded that the value ai/a^ = 1.10-1.15 observed in Baade 
window (1, b) ~ (1°,— 4°), with which our measurements 
are consistent, can be explained by rigid rotation. The pres- 
ence of any disk stars, however, will also contribute rota- 
tional broadening to ai. Since in the vicinity of our fields 
the disk fraction increases closer to the plane, it follows 
that the measured gradient from equation Q could be due 
to the disk contamination. The changes of skewness in the 
/ii distribution tend to support this (see Fig. jHJ. Another 
possibility is that the rotation rate of the bulge a c tually 
increases at lower [6|, as found bv llzumiura et alJ dlQOSh 
from the radial velocities of 124 SiO masers in the Galactic 
bulge. It has been observed that for giants in Baade win- 
dow the metal-poor stars display more spread in the verti- 
cal motio n and less anisotropy when compared to metal-rich 
samples llZhao. Soergel fc RichI 1 19941: IZhao. Rich fc Bielld 
1996). Both of these metallicity dependencies are quite 
steep, so it is likely that the gradient from equation is 
related to a changing mix of populations with more metal- 
poor stars closer to the Galactic bulge minor axis. 

We are not aware of any previous detections of the 
cross ter ms in the Galactic bulge velo city field except the 
report bv lZhao. Soergel fc RichI il994l) of a significant ver- 
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tex deviation between the radial and longitudinal motions 
from Cri. That result is based on a photographic sampl e 
of ~200 K and M giants from ISpaenhaxier et al.1 l|1992ll . 
We note that the latter sample actually shows a hint of a 
sli ghtly negative covariance be tween ^i and ^b (c.f. Fig. 1 
of IZhao. Spergel fc RichI 119941) . The superb resolution of 
HST enabled very significant detections of the Cib cross 
term in many fields. The non-diagonal elements of the 
velocity tensor are crucial to determining the dominant 
orbit families, the importance of streaming motions and 
the need for the intrinsic anisotropy versus solid bo dy ro- 
tation in the Galactic b ulge (Zhao, Sporecl & Ric;h|_^2a 
IZhao. Rich fc Biell o" 1996; Hafner et al.ll2000t l. lHafner et alJ 
(|2000|) published detailed calculations of Cm for several 
lines of sight at positive longitudes including Baade window 
(l°,-4°,Cib = 0.04), and two other: (8.4°, -6°, Cib = 0.15) 
and (1.21°, -1.67°, Cib = 0.04). Taken at face value these 
predictions are roughly of the same magnitude as the results 
from §|1| but have the opposite sign. For a proper comparison 
wit h dynamical models lik e the ones in lHafner et alj (1200(1) 



and lSissantz et al.l (|2004fl we need to wait until the calcu- 
lations are folded with the appropriate selection functions, 
since our measurements are based on substantially deeper 
data than most of the previous samples. 



6 SUMMARY AND CONCLUSIONS 

The main results of our proper motion mini-survey are: (1) 
high quality proper motion measurements for hundreds of 
stars in 35 lines of sight across the Galactic bar, (2) estab- 
lishing the presence of spatial gradients in dispersions cri, CTb 
and the amount of anisotropy ai/ah, and (3) the first reli- 
able detection of the covariance term Cib of the transverse 
velocity tensor. We cross-validated our measurements with 
the ground based O GLE-II data of ISumi et alJ (|20o3) and 
a benchmark study of lKuiiken fc RichI (12002^ . The observed 
slow rise of cti toward the Galactic plane is likely due to the 
increasing disk contamination and/or a possible gradient in 
the bulge rotation speed. The increase in a"b toward the mi- 
nor axis of the bulge is accompanied by the decreasing ratio 
CTi/cTb and possibly results from increasing fraction of metal- 
poor stars. Another possibility is that ah increases due to the 
larger surface density of stars at low I (closer to the Galactic 
centre). We clearly detect the covariance term Cib ~ —0.10 
that implies a significant tilt of the Galactic bulge veloc- 
ity ellipsoid with resp ect to the Galactic p l ane. Using th e 
same procedures as in lBinnev fc Merrifieldl il998. S10.3.2f) . 
we find the tilt is roughly equal to —24°. 

The data presented in this paper provide qualitatively 
new constraints on dynamical models of the inner Galaxy 
and dramatically improved number statistics. Furthermore, 
it may be possible in the near future to augment our proper 
motion samp les with the distance an d metallicity estimates. 
As shown bv lKuiiken fc RichI (120021) . deep color-magnitude 
diagrams can supply sufficiently accurate distance informa- 
tion to effectively isolate the bulge population. In order to 
maximize the discriminating power of model comparisons 
the focus should be on extending the coverage to negative 
longitudes and locations further from the Galactic centre. 
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